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Abstract 

We study unitary invariant random matrix ensembles with singular potentials. 
We obtain asymptotics for the partition functions associated to the Laguerre and 
Gaussian Unitary Ensembles perturbed with a pole of order k at the origin, in the 
double scaling limit where the size of the matrices grows, and at the same time 
the strength of the pole decreases at an appropriate speed. In addition, we obtain 
double scaling asymptotics of the correlation kernel for a general class of ensembles 
of positive-definite Hermitian matrices perturbed with a pole. Our results are 
described in terms of a hierarchy of higher order analogues to the Painleve III 
equation, which reduces to the Painleve III equation itself when the pole is simple. 


1 Introduction and statement of results 

We study unitary invariant random matrix ensembles on the space T-L^ of n x n positive- 
definite Hermitian matrices defined by the probability measure 

— (detM)“ exp [—nTr Ffc(M)] dM, a > —1, (1-1) 

with 


n 

dM = rj (JRe dRe MijdlmMij, (1-2) 

j=l 


and 

Cn= f (detM)"exp[-nTrVfc(M)]dM, (1.3) 

Jut 

in the case where the potential 14 (x) has a pole of order k, i.e. 

Ffc(x) = F(x) + . (1.4) 

The regular part V of the potential is a real analytic function on [0, -|-oo) subject to 
some constraints which we will detail later. In particular, we will require V to be such 
that, for t = 0, the limiting mean density of the eigenvalues as re —>■ oo is supported on 
an interval of the form [0, b], with 6 > 0. 
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It is well known that the eigenvalues of a random matrix drawn from the ensemble 
ra form a determinantal point process. The joint probability distribution of the 
eigenvalues is given by 

n 

— —A(x)^ A(x) = {xj-Xi), (1.5) 

j = l 

with the partition function Zn^k = Zn^kiY) given by 


Zn,k= / ( 1 . 6 ) 

./[0,+oo)- 

For t > 0 small, the model (|l.ip can be seen as a singular perturbation of the unitary 
invariant ensemble corresponding to t = 0. For t > 0, the eigenvalues are pushed away 
from 0 because of the pole in the potential, and for large n, the probability of finding 
eigenvalues close to 0 is small if t > 0 is independent of n. However, if t —>• 0 together 
with n —)■ oo, this repulsion becomes weaker and one expects a transition between a 
regime where eigenvalues are likely to be found in the vicinity of the origin, and one 
where eigenvalues are unlikely to be found near the origin. 

The effect of singular perturbations of unitary invariant ensembles has been of recent 
interest [3 H EZl EH]. In m, the singularly perturbed Laguerre Unitary Ensemble 
(pLUE) was studied, given by the measure 


^(detM)“exp 


—n Tr 



dM, 


(1.7) 


on the space T-L^i where Cn is a normalisation constant. A relation between this model 
and the Painleve III (henceforth PHI) equation was established for fixed n in [5]. In 
subsequent work [4|, a singular perturbation of the Gaussian Unitary Ensemble, which 
we will refer to as pGUE, was studied, defined by the measure 


1 

exp 

a. 


, 1 n t 

-n Tr -M^ H- 

' 2 2M2 


dM, 


( 1 . 8 ) 


on the set of n x n Hermitian matrices T-Ln- 

In [3] the double scaling limit where f —>• 0 as n —)• oo of the partition function was 
analysed using Riemann-Hilbert (RH) techniques. A connection to PHI was also found 
by relating the orthogonal polynomials associated to the pGUE measure to those of the 
pLUE measure. In [27], the double scaling limit for the eigenvalue correlation kernel 
in the pLUE was studied. A limiting kernel was found, defined in terms of a model 
RH problem associated to a special solution of the PHI equation. This limiting kernel 
degenerates to the Airy kernel if t approaches 0 at a slow rate, and to the Bessel hard 
edge kernel if f —)• 0 at a fast rate. In [28], asymptotics for the partition function in the 
pLUE, again in terms of a PHI transcendent, were obtained. 

In the present work, we will obtain asymptotics for the eigenvalue correlation kernel 
for a fairly general class of potentials V, perturbed with a pole of order A: G N, in 
a double scaling limit where the strength of the perturbation, t, goes to zero at an 
appropriate speed as the size of the matrix, re, is taken to infinity. The double scaling 
limit will be tuned in such a way that the macroscopic behaviour of the eigenvalues in 
the large re limit is the same as in the non-singular case where t = 0, but such that 
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the microscopic behaviour of the eigenvalues near 0 is affected by the singularity of 
the potential. In addition, we will obtain double scaling asymptotics for the partition 
functions associated to the LUE perturbed with a pole of order k, and to the GUE 
perturbed with a pole of order 2k. Our results will be described in terms of special 
solutions to a family of systems of ODEs, indexed by A: G N, which can be seen as a 
hierarchy of Painleve III equations. 


1.1 Motivations 

There are a number of motivations for considering the model (II.II) . Firstly, it was ob¬ 
served that critical one-matrix models in which the limiting mean density of eigenvalues 
vanishes in the bulk of the spectrum or at a higher than generic order at the edge of 
the spectrum, are in one-to-one correspondence with (super-)Liouville field theories. In 
this context, it is natural to look for models with higher order critical points in which 
multiple zeros appear in the bulk or at the edge of the spectral density, as it is known 
that such higher order critical models correspond to coupling a minimal conformal field 
theory to (super)-Liouville theory. A review of these facts in the non-supersymmetric 
case can be found in [12]. Supersymmetric versions followed later in [20| and are nicely 
reviewed in the appendix of [24]. One-matrix models in which the potential has poles 
exhibit a new type of critical behaviour. Although at present no conformal field theory 
analogue of such models is known, it is again natural to study higher order critical 
points in this context. 

A second physical motivation arises in the field of quantum transport and electrical 
characteristics of chaotic cavities. Here the quantity of interest is the Wigner-Smith 
time-delay matrix Q, the eigenvalues of which, tj, are known as the “partial delay 
times”. In systems in which the dynamics is chaotic, a RMT approach has been quite 
successful and one of the central results of this approach is the joint probability density 
for the inverse delay times, jj = first obtained in [BJlT]. It takes the form 

P(7i,...,7„) = |A(7)|^]j7|"/^e"2'^R (1.9) 

j=i 

where /3 depends on the symmetries of the system, with /3 = 2 a common case. Since 
many observables may be expressed in terms of Q, the problem of computing expecta¬ 
tion values with respect to the above measure is relevant. In particular the observable 


I 1 

Tw = -TtQ = - 


n 


n 


i=l 


(I.IO) 


known as the Wigner time-delay has been considered in the recent work [25]. The 
partition function for our model (II.II) coincides with the moment generating function 
for the probability density of tw in the case A: = 1 and a = /3n/2. Although we will not 
scale a with n in the present paper, the moment generating function in this model is the 
partition function for a matrix model which shares the feature of a singular potential 
with our model (|I.1I) . 

Furthermore, in a very recent paper [18], an observable Rq known as the “charge 
relaxation resistance” defined as 


Rq oc 


TrQ2 

(TrQ)2 




( 1 . 11 ) 
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was considered. To compute Rg, the approach taken in [TH] was to compute the distri¬ 
butions of TrQ^ and TrQ separately. The associated moment generating function is 
the partition function for a perturbed pLUE model with a singularity in the potential 
of order k = 2, thereby demonstrating the physical relevance of the model studied here 
for A: > 1. 

A third motivating model appears in the field of spin-glasses [2]. Here a model 
corresponding to the GOE perturbed by a pole of order k was analysed for its relation 
to the distribution of the spin glass susceptibility in the Sherrington-Kirkpatrick (SK) 
mean-field model. The GUE version of the partition function in such a model relates 
directly to the partition function in the pLUE, as we will see later on in this paper. 
One-matrix models in which the potential has poles have also appeared in the context 
of replica field theories [21j . 

Einally, our last motivation for this work is to seek a natural candidate for a PHI 
hierarchy. The notion of a PHI hierarchy has appeared little in the literature; one of the 
few mentions being [23]. This work was partly motivated by the desire to understand 
whether the hierarchy proposed in [23] would appear when the order of the pole was 
increased and if not, what form the alternative hierarchy takes. It appears that the 
hierarchy of equations which we will obtain is different from the one in [23] . 

1.2 Statement of results 

Our main results are the following. 

1. We define a hierarchy of higher order PHI equations and prove the existence of 
special pole-free solutions to it. 

2. We obtain double scaling asymptotics for the partition function in the LUE per¬ 
turbed with a pole of order A: > 1, in terms of a higher order PHI transcendent. 
This generalizes the result from m for A; = 1. 

3. We obtain double scaling asymptotics for the partition function in the GUE per¬ 
turbed with a pole of order 2k, k > 1. They are again given in terms of higher 
order PHI transcendents. In the case A: = 1, such asymptotics were already 
obtained in [3], but written in a different form. 

4. In the double scaling limit, we prove, for A; > 1 and for general V, that the 
eigenvalue correlation kernel near 0 of the model (II.ip tends to a limiting kernel 
built out of a model RH problem associated to the PHI hierarchy. This extends 
the result from m for V{x) = X and A: = 1. 

We now state our results in more detail. 

A Painleve III hierarchy 

Given k G N, consider the system of A; -|- 1 ODEs indexed by p = 0,..., /c, 
p 

'y ^ {^k—p+q+l^k—q ~ i^k—p+q^k—q) + ^^k—p+q^k—q ~ ^'^^k—p+q^k—q) = Tp, ( 1 . 12 ) 
q=0 

for A: -|- 1 unknown functions u = u{s) and ii = ii{s),... ,ik = ^k{s), with 

4+i(s) = o, 4(s) = |. (1.13) 
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The Tp’s are real constants that play the role of times. The p = 0 equation always 
results in 

^ ^ ” 4 ^" ~ + ''“o) • (1-14) 

Substituting this expression for u in the other equations, we are left with k equations 
for k unknowns £i,... We refer to this system of equations as the k-th member 
of the Painleve III hierarchy. Eliminating for 2 < p < k leads to an ODE for ii of 
order 2k. 


Example 1.1 Eor fc = 1 we have the equation 

/TO = _ '^0 , D 

^ ll{s) S S ^l{s) s ’ 

which we identify as a special case of the Painleve III equation, see m- 


(1.15) 


Example 1.2 If A: = 2 we have a system of two ODEs; 

n To , i 2 {s)'^ ^l{s)i 2 {s) , ii{s) l'^{s) i 2 {s) 


+ 


+ 


2^i{s)i2{s) i2{sY (■2 {sY h{s)(.2{s) h{s) (■2{s) 2 £i(s) 


= 0 , (1.16) 


and 


h{sf^'2{sY_^,^^^^2j^sl'2{sf Toli{sf STq 


l2{s? ' Us 

2Us)H'i{s) 


-m- 


T2 


Us) 


Us)"^ Us) 

- 2Us)Us) + sUs) + 2Us)Us)- 


(1.17) 


One can eliminate £2 in order to obtain a single equation of order four for ^ 1 . 


We can construct, for any fc = 1, 2,..., a special set of solutions to the 

k-th member of the PHI hierarchy in terms of a model RH problem. The function l'i(s) 
will be of particular importance to us. The model RH problem consists of finding a 
function T = T(z; s) satisfying the following properties. 


RH problem for T 

(a) T : C \ S —analytic, with S = U {0} as illustrated in EigurelH The 

half-lines Si, S 3 can be chosen freely in the upper and lower half plane. 

(b) T has continuous boundary values 4/^(2:) as 2 : G S \ {0} is approached from the 
left (-I-) or right (—) side of S \ {0}, and they are related by 


T+( 2 ) = T_( 2 ) 
T+(z) = T_(z) 
T+(z) = T_(z) 


1 0 

Trio. 


—e 

0 -1 
1 0 


—e 


1 0 
—Tria 2 ^ 


z G Si, 
2 G S2, 

2 G S 3 . 


(1.18) 

(1.19) 

( 1 . 20 ) 
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Figure 1: The jump contour F for the model RH problem for T. Contours are labelled 
by S and sectors by hi. 


(c) As z —>■ oo, there exist functions p{s),q{s),r{s) such that 'h has the asymptotic 
behaviour 


1 / qis) -ir{s) 
z \ip{s) -q{s) 




( 1 . 21 ) 


where ) > with ui ~ where the principal branches of 

and z~^l^ are taken, analytic off (—oo,0] and positive for z > 0. The third 
Pauli matrix us is given by <73 = 


0 -1 

(d) As z —)• 0, there exists a matrix 'I'o(s), independent of z, such that T has the 
asymptotic behaviour 


T(z) = To(s)(/ + O(^))e"("-) 

( 1 . 22 ) 

G Hj, where Hi,H 2 , H 3 are given by 


II 

(1.23) 

^2 (_gTO 1 ) . 

(1.24) 


(1.25) 


Remark 1.3 In the case k = 1, the RH problem for T coincides, up to the orientation 
of the contours, with the model RH problem introduced in [27j . In the case A; = 0, it is 
a RH problem which can be solved using Bessel functions [26]. 


6 



Remark 1.4 It is important to note that the function 




(1.26) 


is not analytic at z = 0 and in fact has a jump across r 2 . Indeed, if the asymptotic 
behaviour in condition (d) of the above RH problem were modihed to 

(1.27) 


'I'(z, s) ='I'o( 2 :, s)e ^ 


with ^o(2) s) analytic at z = 0, then the resulting RH problem would have no solution. 
A more detailed description of the analytic structure of ^ near the origin can be given 
as follows. Define 


ft \ ^ f \ \a u du 

f 2 {z,s) = — / |u|"e“e -1 - 

27 ^^ Jq 


u — z 


and let be defined by 
^{z, s) = ^o{z, s) 




(1.28) 


(1.29) 


Then Tq is an analytic function near z = 0, and Tq dehned by (jl.22p takes the form 


^o(s) = 4-0(0, s) 


1 /2(0,s) 

0 1 


(1.30) 


Remark 1.5 If the model RH problem has a solution, it follows from standard tech¬ 
niques that the solution is unique. Existence of a solution is a much more subtle issue. 
We will show that the model RH problem is solvable for A: £ N and s > 0. 

Theorem 1.6 Let a > —1, and let 'I'(z;s) be the unique solution of the model RH 
problem for s > 0. Then, the limit 


Vais) =—2i-^ lim z'I'(z, s^)e ^z*^^ 

ds U->-oo 




(1.31) 


exists and it is a solution of the equation for t\ in the k-th member of the Painleve III 
hierarchy, with the parameters tq, ... ,Tk given by 

42fc-|-1^2, — 

Tp = < forp = k, (1-32) 

0, for 0 < p < k. 

Moreover, y{s) = yais) has the following asymptoties as s ^ -|-oo and as s —)• 0, 


vis) = 


8k 


3 

2 


2k+ 1 

2/(s) = 0(s2^-1) + 0(s2“+^), 
where we have used the constants 
2^-^{k - 


zo = 


{2k-l)\\ 


ft = 

with (2j + 1)!! = {2j + l){2j — 1)... 3.1 the double factorial. 


+ o{i), 

S +00, 

(1.33) 


s —)* 0, 5^0, 

(1.34) 



(1.35) 



(1.36) 
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Remark 1.7 In fact we will prove a more general result in the sense that all ij’s may 
be extracted from the model problem for To this end we write the asymptotic 
expansion of 'h as 2 ; —)• oo as 


^'( 2 ) = 1^1 +1 

Define the formal power series 


^ W.(«) 


m2{z,s) = 


i=i 


rj{s^) 


(1.38) 


and dehne in addition formal power series mi, m 3 , mi in z mi{z, s) = ^■ 

The quantities m 3 and mi are defined in terms of mi and m 2 as 

r (s^) (mi + 1 ) 1 

2 * 


m‘i{z,s) := —dgmi - —- — -8^1712 


/ 2^ ^ ' 3r (s2) \ 

+ ’"2(3r - + (1 


1 fr{s^) 


39) 

mi( 2 , s) := ^ {^^^^1712 - ds7n2'j , (1.40) 

and mi can be found recursively from the relation 
1 1 

= 2 ("^4 + " 12 ^ 3 )^ - 2 X] (1.41) 

i=0 

where (m| + 17127173 )j denotes the coefficient of z~^ in the formal power series of m| + 
17121713 . Note that upon substituting in the expressions for m 3 and mi in terms of mi 
and m 2 , the right hand side of (|1.41l) only contains rrii^i with i < j and therefore gives 
a well defined recursive relation for mij. 

We now introduce the matrix 


K{z,s) := 


mi + mi + 1 m 2 

m3 + s“^r(s^)(mi + mi + 1) mi — mi + s“^r(s^)m2 + 1 

where r = ri. We then obtain expressions for all ij in terms of rj via the relation 


, (1.42) 


4 


zf = 


(1.43) 


n=0 


(4z)^+^ Tr 


1 


{dzK)K-^ + -Ka-K-^ + —K ( scj+ - -CJ 3 ) 


1 


(1.44) 

where the notation (...)+ denotes a projection onto the positive power parts of the 
power series. The above relation gives for the first few ^j’s, 
s 


4(s) = 


.^i(s) = —4sr' (s^) 
8 
s 


(1.45) 

(1.46) 


■^2(5) = - {r (s^) (s^) + 2 s^r' (s^) ^ — 2 r (s^) (/ (s^) + s^r” (s^)) — 2r2 (s^)) • 

(1.47) 


Note that these expressions are independent of k. 











Remark 1.8 There has appeared a distinct definition of a Painleve III hierarchy in 
the literature [23]. The model problem for ^ may be connected with this alternative 
hierarchy if we choose a different time to act as the independent variable in the ODE. 
To see this, we generalise the model problem by altering the behaviour at the origin 
to 

^{z) = (1.48) 

and then define 

(1.49) 

Using the asymptotic behaviour of ^ in the definition of the Lax pair 

and 1 := (1.50) 


we find 


B = 





(1.51) 


and A = where Aj are si dependent matrices and u and v are some 

undetermined functions. The form of these Lax matrices precisely matches those pro¬ 
posed in [23j as forming a Lax pair for a Painleve III hierarchy. However, in the random 
matrix model under consideration, it is more natural to work with the variable s instead 
of si, and this leads us to a different Painleve hierarchy. 


Double scaling limit for the partition functions in the perturbed LUE and 
GUE 


Define the pLUE partition function as 




■pLUE(,~, _ 


n,k,a 


{t) = 


1 


AixfH 


./[0,-|-cxd)'‘ 

and the pGUE partition function as 


xje 


—n[ Xj+[ — 


dxj, a > —1, 






(1.52) 


dxj, a > -1/2. (1.53) 


Whereas for the pLUE and pGUE partition functions with k = 1, large n asymptotics 
are described in terms of a special solution to the PHI equation [U [28| , when perturbing 
a unitary invariant ensemble with a pole of order /c > 1, the special solutions to the 
higher order analogues of PHI introduced before will appear. 


-- ^*^+1 

Theorem 1.9 Hs re — )• oo and t — )• 0 m such a way that Sn,t ■= “2 fc t ^ s, we 
have the asymptotics 




n,k,0' 


2 r^n,t 




X (l + 0(re-^)), (1.55) 
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where r^is) := r(s; a) is related to the Painleve transcendent ya{s) by 


ya{s) = -2^ (r’a(s^)) , r„(0) = ^(1 - 40^). 

In other words, we have 

1 9 1 

^ais) = -(1 - 4a2) - 2 ya{v)dv. 


(1.56) 


Double scaling limit of the correlation kernel 

In what follows, we consider potentials V which are real analytic on [0, +oo) and such 
that 


Vix) 


lim 

x^+cxD log(a:^ + 1) 


= + 00 . 


(1.57) 


The correlation kernel for the eigenvalues in the model o is expressed using 
orthogonal polynomials with respect to the weight 


w{x) = x"e = x“exp 


—n 


^l"(x) + 



(1.58) 


on [0, +oo). Let pj, j = 0,1,... he the family of monic polynomials of degree j charac¬ 
terised by the relations 

Pj{x)pm{x)w{x)dx = hjdjm- (1.59) 

The correlation kernel of the determinantal point process (11.51) can be written as 



Kn{x,y) = (p^{x)pn-l{y) - Pniy)Pn-lix)) . (1.60) 

X y 

Important quantities that can be computed directly from Kn{x,y) are the /c-point 
correlation functions. The one-point function, or limiting mean eigenvalue density, is 
given by 

p{x) := lim —Kn{x,x), (1-61) 

n^oo Tl 


and describes the macroscopic behaviour of the eigenvalues in the large n limit. For 
fixed t > 0, the measure du^x) = p{x)dx is characterised as the equilibrium measure 
which minimizes 




- -dv{x)dv{y) + 

x-y\ 


Vk{y)dv{y), 


(1.62) 


among all Borel probability measures v on [0,-|-oo). If we let n —?> oo and at the 
same time t —?■ 0, the limiting mean eigenvalue density is not affected by the singular 
perturbation and is equal to the equilibrium measure dp = '4){x)dx obtained from 
minimising 



- -dp{x)dp{y) + 

x-y\ 


V{y)dp{y), 


(1.63) 
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in which the singular part of 14 lias been dropped. It is this measure and density we 
will use in the remainder of the paper and refer to as the equilibrium measure. The 
measure // is characterised by the Euler-Lagrange equations [22], stating that there 
exists ^ G M such that 


2 

2 


log \x — y\’ip{y)dy — V(x) = i, x £ suppy, 
log \x — y\'tp{y)dy — V(x) < £, x £ [0, +oo). 


(1.64) 

(1.65) 


We will require that V is such that the equilibrium measure y is supported on a 
single interval of the form [0, 6], which implies [8| that its density 4 can be written as 


ijj{x) 





X £ [0,6], 


( 1 . 66 ) 


where /i is a real analytic function, non-negative on [0,6]. We will also require that V 
is regular, in the sense that 


(i) the density 4 is positive in the interior of its support (0,6), i.e. h{x) > 0 for 
0 < a; < 6, 


(ii) 4 vanishes like a square root at the endpoint 6 > 0 and like an inverse square 
root at the endpoint 0, i.e. 6,(6), 6(0) / 0, 

(iii) the variational inequality (|1.65p is strict on (6,+oo). 

Define functions 'ipj(x',s) in terms of the solution to the model RH problem for 'k 
as follows. For £ Di, let 


\i/j2{z,s)) ■ 


4'(z,s) 



(1.67) 


and for z G D2 U D3, we define and "02 in such a way that they are analytic in C \ Si. 
This means in particular that 


\'tjj2{x,s)J 


4'+(x,s) 


—e 


1 

—'Kid 


( 1 . 68 ) 


for a; < 0. 

Then, ipi,'4’2 satisfy the system of equations 


dz 



if d — ^bs ■2U{s) 

s \is^{c + ds~^U{s) — ^bs~^U{s)‘^) 


— IS 


— CL + 



(1.69) 


and they are characterised as the unique solution with the asymptotic behaviour 





(1.70) 


as z —)• 00 in (and, more generally, for j arg^j < vr — e for any e > 0), and 



(1.71) 
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1 

as z —)• 0 for 2 ; e Here, U{s) = y{s')ds', d = a(s“^ 2 :, sz), b = 6 (s“^z,s 2 ), 

c = c(s~^z, sz), with 


a{z,s) = --dsb{z,s), 

c{z, s) = {z- u)b{z, s) - ^d^b{z, s), 

and 


(1.72) 

(1.73) 


b{z,s) = ^ ^4-n(s)(4^) 

^ ' n =0 


(1.74) 


Theorem 1.10 In the double scaling limit where n —)> oo and simultaneously t —)• 0 in 


1 2fc+l 


such a way that 2 fccin k t ^ s > 0, with ci = bh^Of, we have 


1 / -u -V 

hm - 

n^oo c\n^ V cin^ cin^ 


t = 


^Piii 


{u,V]s), 


for u,v < 0, where 


V] s) = e!' 


^'ilji{u;s)ij 2 iv;s) - ijiiv, s)'ilj 2 iu-, s) 


(1.75) 


(1.76) 


2'Ki{u — v) 

The limit is uniform for s in compact subsets of (0, +oo). 

Remark 1.11 For k = 1, our limiting kernel {—u, —v; s) is equal to the limiting 

kernel obtained in m- 

Theorem 1.12 The limiting kernel {u,v] s) is positive for u,v < 0 and s > 0, 

and it has the following limits, 

(1.77) 

.•s’’= A(n,n), (1-78) 

where rj = Sa{u,v) and A{u,v) are the Bessel and Airy kernels defined as 


lim v; s) 

s^O 

= Jq(- 

-u, 

-v) 

i, 





Zt, 

/ / 




lini 

l^zo 

+ 

V 

S 3 

u 

s—>-+oo C2 



C2 


Iq.{u,v) = 
A{u, v) = 


Ja{Vu)y/vJa{^/v) - Ja{y/v)y/uJ’^{y/u) 
2{u — v) ' 

Ai (n)Ai'(?;) — Ai (n)Ai'(n) 


u — V 

The constant C 2 has the explicit form 


(1.79) 

(1.80) 


C2 = 




j=0 


2Jj! 


(1.81) 
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Outline 


In Section [21 we will study in more detail the model RH problem for 'h. Using Lax pair 
type arguments, we will establish the relation between the RH problem and the PHI 
hierarchy. In addition, we will prove the solvability of the RH problem for s > 0, which 
implies that the PHI solution y{s) has no singularities for s > 0. In Section [3l we will 
study the large n asymptotics for the orthogonal polynomials associated to our random 
matrix model by means of the Deift/Zhou steepest descent method. A crucial feature 
in this analysis will be the construction of a local parametrix near the origin in terms 
of the model RH problem for In Section HI we will obtain large s asymptotics for 
'I' and for the PHI solution y using an asymptotic analysis of the model RH problem. 
In Section [5l small s asymptotics for 'h and y are obtained. The results from Section 
[2] together with the asymptotics from Section HI and Section [5| will lead to a proof 
of Theorem [RHI In Section [6l we use the asymptotics for the orthogonal polynomials 
obtained in Section [3] to prove Theorem ll.lOl In Section [71 we use a differential identity 
which enables us to derive asymptotics for the pLUE and pGUE partition functions 
from the asymptotics for the orthogonal polynomials, and to prove Theorem 11.91 


2 Properties of the model RH problem 

2.1 Connection to the Painleve III hierarchy: proof of Theorem 11.61 

We begin by introducing a new model RH problem for a function $, obtained from the 
model problem for defined in the introduction. Let 

( 2 . 1 ) 

with r defined by ()1.2ip . Using the RH conditions for T, we find that satisfies the 
RH problem given below. 


RH problem for <I> 

(a) <I> : C \ S —^ analytic. 

(b) <1 has the jump relations 


^+(2) = ^-(^) ( Ji, 


$+(z) = $_(2) 


0 -H 
1 0 / ’ 


^+(^) = ^-(2) ( Jyi 


(c) As 2: —)• 00, <1 has the asymptotic behaviour 

1 / w{s) v{s) 


<^{z) = 


1 0 
v{s) 1 


1 + 


Z \h{s) -w{s) 


z G Si, 

(2.2) 

2 G S 2 , 

(2.3) 

z G S 3 . 

(2.4) 

2 - 2 )) 

■1710-3 ^-i(T3 ^ 

(2.5) 
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where N = + iai) as before, and 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


v(s) = s ^r(s^) 


h{s) = s ^p(s^). 


(d) As z —)• 0, there exists a matrix ‘ho('S), independent of z, such that has the 
asymptotic behaviour 



^z) = Ms){I + 0{z))e 


(2.9) 


for z in sector Qj, where Hi, H 2 , are given in the model problem for 'h, see 




Under the assumption that the model RH problem is solvable for s > 0 (we will prove 
this assumption in Section [2.2jl . we are now in a position to prove Theorem 11.61 First 
we note that the limit in (ll.3ip exists by a simple application of (ll.2ip . and that the 
differentiability with respect to the parameter s can be proved using general methods 
which involve re-writing the RH conditions as a singular integral equation. We have 


y{s) = - 2 ^(r(s 2 )) = -2-^{sv{s)). 


( 2 . 10 ) 


For the relation to the PHI hierarchy we use well-known methods of isomonodromic 
deformation theory. In particular we define the Lax matrices, 


A := (d^<P) and B := (d,^) 


( 2 . 11 ) 


and note that, due to the jump matrices of $ being constant, the Lax matrices are mero- 
morphic functions in z, and may only have isolated singularities. Using the asymptotic 
expansion of near z = 00 and z = 0 one can easily show that A and B are Laurent 
polynomials of the form 


fc-i-i 




( 2 . 12 ) 


j=0 



(2.13) 



Computing the (l,2)-entry of the 


z ^ term in B using the asymptotic behaviour of <I> at 2 ; = 00 together with the fact 
that B is regular at zero implies 


v'{s) — v(s)'^ + 2w{s) = 0 . 


(2.14) 


We can also give explicit expressions for the elements of Aj in terms of coefficients 
appearing in the asymptotic expansion of near 2 ; = 00 . Although this will prove 
largely unnecessary, one important result obtained in this way is that 


(^ 1)12 = I 
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(2.15) 

(2.16) 











where in the second line we have used (j2.14l) . To complete the proof we now show 
that u{s) is the function appearing in ()1.12l) . To this end, consider the compatibility 
condition, 


dzds^ = dsd^^, 


(2.17) 


which, upon substituting in the expression (12.111) for the derivatives in terms of the 
Lax matrices, becomes. 


d,B-dsA+[B,A] = 0. 
If we parameterise A as 


A{z,s) = 


a{z,s) b(z,s) 
c{z, s) —a(z, s) J ’ 


then (I2.18P is equivalent to three coupled ODEs, 


a{z,s) = --dsb{z,s), 

c{z, s) = {z- u)b{z, s) - ^d^b{z, s), 
dsc{z, s) = 1 + 2{z — u{s))a{z, s). 


(2.18) 

(2.19) 

( 2 . 20 ) 

( 2 . 21 ) 

( 2 . 22 ) 


The first two equations provide a and c in terms of b. Although the hnal equation can 
be used to obtain an equation for b, this is not the course we will follow here, however 
we will remark upon this possibility later. 

So far we have only used the asymptotic behaviour of ^ near zero to determine the 
order of the pole occurring at the origin in each Lax matrix. To complete the proof we 
now consider the behaviour of <1> near zero in more detail. Working in Hi, see Figure 
m write (lOl) as 


$(2) = $o(s)i^(2)e^(^^‘"^ 

where F{z) = I + 0{z) as z ^ 0, and A can be written in the general form 


A(2:) := -Aologz + ^ 


A, 


n=l 


nz'‘ 


Ao = Afc = Ai = ... = Afc_i = 0 . 

Note that det A can be written as 


(2.23) 

(2.24) 

(2.25) 


det A = det (S^^xh ^) = — A'(z)^det ^1 + ^ a^F'{z)F{z) , (2.26) 


and that as 2 ^ 0, 

A'(2) = 0(z-"-1), 

F'{z)F{z)-^ = 0(1). 

Hence we obtain the estimate 

detA = -A'(z )2 + 0(2-^-1), 


(2.27) 

(2.28) 


(2.29) 
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as 2; ^ 0. Using (|2.2UI) and (12.211) and rearranging gives 


+ (2 - K(»))(>" - - 


A. 

dz 


A(z)) =0{z->^-^), 


0 , 


(2.30) 


where primes denote derivatives with respect to s. The important observation is that 
the left hand side behaves as 0{z~‘^^~‘^) as z —>■ 0 which after comparing with the right 
hand side yields k + 1 equations. Explicitly, in (I2.30j) . let 


b{z,s) = 


(42)fe+i 


^4-i(s)(4^ 


with Iq = s/2. This gives 

k 


(4_2,)2fc+2 


E 

j,m=0 


[{ik-A-mf - K-/k-m + 4u4-/fc-m] 
= + 0{z ^ ^), z —)• 0. 


(2.31) 


(2.32) 


After some rearrangement and dropping terms of order 0{z ^ ^) the above equation 
can be written as 

4 ^ P 

(A '\2k+2 ^ ^ y Ak—p+q+l^k—q ~ i^k—p+q^k—q) + ^^k—p+q^k—q ~ ^'^^k—p+q^k—q) (4z)^ 

p=0 q=0 

= (^^A(z)y+ 0(z-'=-i), (2.33) 

Finally, substituting in (I2.24p and comparing coefficients yields (|1.12l) . where Tp is a 
constant taking the form 

p 

Tp = A-P+^ Y, A-p+qAk-q. (2.34) 

q=0 

In our case, by (|2.25l) . we have (ll.32p . This shows that u{s) appearing in (12.131) together 
with the quantities Ij satisfy the system of equations (jl.l2p . The relation between 
li{s) and v{s) can be obtained from (I2.10p and (I2.16P together with the fact that 
(^2)12 = 61 = \li. 

We now prove the statements appearing in Remark 11.71 We parameterise the be¬ 
haviour of 4* as z ^ 00 by 


4>(z) 



0 \ (j , (mi + 1714 

V V V ^3 





(2.35) 


where m* are formal power series in z ^ and we note that this can be written as 


4>(z) = 


(2.36) 
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Computing the Lax matrix B using the above expression yields the relation 

z[a-,K] = dsK - 2v'{s)a-K - [a+,K], (2.37) 

which upon substituting in v{s) = r{s^)/s gives ()1.39p and (ll.40p . Furthermore the 
fact that detiL = 1 yields the recursion relation (ll.4ip . Finally, computing the Lax 
matrix A using (I2.36p gives 

A = (d,K)K-^ + ^Ka.K-^ + (^sa+ - K-\ (2.38) 

from which the Ij may be extracted by projecting out the 1, 2 element of A. 

Remark 2.1 Substituting (I2.20p and ()2.21l) into (j2.22p yields 

zdsb{z, 'S) = ^ s) + Au{s)dsb{z, s) + 2u'{s)b{z, s)) + i. (2.39) 

We may compute b{z, s) by substituting (j2.31l) into the above equation, which yields 

= q + 4u^'. + 2n'£j, (2.40) 

for j > 0 subject to the conditions = s/2 and = 0. By iterating the above 
equation one obtains each Ij as differential-integro polynomials in u. The final condition 
^k+l = 0 yields an ODE for u of order 3A; +1, which is an equivalent way to characterise 
the A:-th PHI equation. The reason we prefer ()1.12l) is that it yields an ODE for ii of 
order 2k which coincides with the third Painleve equation for k = 1. 

Einally, it is interesting to note that (|2.40l) is precisely the recursion equation for 
the Lenard differential polynomials appearing in the definition of the Painleve I and 
II hierarchies. The difference between the Painleve I, II, and III hierarchies lies in the 
initial condition; in the PI case it is Iq = —4tt, for PII = 1/2, and for PHI = s/2. 
This difference in initial condition has a large effect — in the case of PI and PII the 
left-hand-side of (j2.4np is always a total derivative and therefore the recursion relation 
can be integrated at each order; this is not the case for the PHI hierarchy. 

Recall the definition of 'i/i, ■02 in (Il-h7|) . 


Corollary 2.2 The ij;-functions associated to T, defined in (ll.67p . satisfy the ODE 
system 


dz 



1 / a — ^bs 2 [/(s) —is 2 b \ 

s \is 2 (c-\-as~ 2 U{s) — jbs~^U{s)‘^) —a-\-^bs~ 2 U{s)j Yl’ 2 ) 


(2.41) 


1 


where U{s) = y{s')ds', a = a{s ^z,s 2 ),b = b{s ^z,s 2 ) and c = c{s ^2:, sz), and 
are characterised as the unique pair of solutions to ()2.4ip with the asymptotic behaviour 

KTOh - Km . 


Proof. Substituting (12.11) in the first equation of (12.lip , and using (12.191) . we obtain 
()2.41l) . The asymptotic conditions are a consequence of conditions (c)-(d) of the RH 
problem for T, and the fact that these determine the solution uniquely follows from 
the general theory for linear systems of ODEs. □ 
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2.2 Proof of existence of ^ 

In order to prove the solvability of the RH problem for 'h, we consider its homogeneous 
version. This consists of the RH conditions (a), (b), and (d), and has condition (c) 
replaced by 

T(2) = 0{z-^) X as z ^ oo. (2.42) 

Lemma 2.3 (vanishing lemma) Let Tq satisfy the RH conditions (a), (b), and (d) 
in the RH problem for T, and the homogeneous asymptotic condition (I2.42p . Then, 
^o(^) =0. 

Proof. This was proved in for k = 1. The proof from [27] is easily extended to 
the case fc £ N: equations (2.20) and (2.24) in |27| have to be modified in the natural 
way, but with the function H{z) as dehned in equation (2.26) in [27|, it is still true that 
H{z) = 0(1) as z —)• 0. The remaining part of the proof relies on this fact, but not on 
a more precise description of the behaviour near 0, and the value of k is unimportant 
here. □ 


For a general class of RH problems, the vanishing lemma (i.e. the fact that the 
homogeneous RH problem has only the trivial solution) implies the solvability of the 
non-homogeneous RH problem. This follows from the description of RH problems in 
terms of singular integral equations. The existence of a solution of the RH problem is 
then equivalent to the bijectivity of a certain operator. This operator is a Fredholm 
operator of index zero and because of this, bijectivity follows from the fact that the 
kernel of the operator is trivial. Triviality of the kernel is equivalent to the vanishing 
lemma. This procedure works for a fairly general class of RH problems, as long as the 
RH problem is equivalent to a RH problem without singular points. The RH problem 
for T has a singular point at 0, but can be transformed to one without singular points 
using the representation (|1.29p . 

Define 


T(z) 


T(z), for l^l > 1 , 
To( 2;), for l^l < 1. 


(2.43) 


Then, T has no jumps inside the unit circle, but it has a new jump on the unit circle, 
which we choose with clockwise orientation, given by 

TiH^)$+(^) = = (J z^^^H^, (2.44) 


on the part of the unit circle in Llj. The RH problem for T has no singular points, 
and the vanishing lemma procedure thus applies to this RH problem, and to the RH 
problem for T. It follows that the RH problem for T is solvable. For more details on the 
vanishing lemma and how it implies solvability of the non-homogeneous RH problem, 
we refer to 13 da la [mill]. 


3 Asymptotic analysis of the RH problem for orthogonal 
polynomials 

The steepest descent analysis we will perform follows closely the calculations in [26] 
(corresponding to the case A: = 0) and in |27| (corresponding to A: = I) which we will 
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refer to extensively. The main difference is in the extra condition imposed at z = 0. 
This will enforce us to construct a local parametrix near 0 which is not built out of 
Bessel functions, but in terms of our model RH problem associated to higher order PHI 
equations instead. 


3.1 Equilibrium measure preliminaries 

As stated in the introduction, a central role in the steepest descent analysis is played 
by the equilibrium measure dfj,{x) = 'il){x)dx. For our purposes, we only need that the 
measure ^ is characterised by the conditions (ll.64l) - (ll.65p . Under our assumptions on 
V, namely the fact that it is such that the equilibrium measure /r is supported on a 
single interval [0,5] and that it is regular, by (ll.66p . the equilibrium density has the 
form 


V>(x) 


1 

27ri 


/i(x)R+(x)y(o,b](x), 


where 


R{z) = 



(3.1) 


(3.2) 


The right endpoint of the support 5 > 0 is a constant depending on the potential V, 
h{x) is a real analytic function in x and X[o,fe] is the indicator function for the set [0, 5]. 
The principal branch of the square root is taken in R, so that R is analytic in C \ [0, 5] 
and positive for z > b. Furthermore, we have h{x) > 0 for x G [0, 5]. 

In the steepest-descent analysis, the following functions will also prove useful. 
Firstly we define the “^-function” 


aiz) 



x)dfi{x), 


(3.3) 


where the principal branch of the logarithm is taken, meaning g is analytic on C \ 
(—00,5]. Under the one-cut regularity assumption which we imposed on V, the g- 
function has a number of properties we will make use of, in particular 

g+ix) + g-{x) - V{x) - £ = 0 xG[0,5], 

g+{x) + g-{x) — V{x) — £ < 0 xG{b,+oo), 

fb 

g+{x) — g-{x) = 27ri / 'ij;{x)dx. 

J X 

Let us also define 

C{z) =-- [ h{s)R{s)ds for 2 G C \ (—oo, 5], (3-7) 

2 Jb 

where the integration path does not cross the real axis. By (13.6p we have, 

g+ix) - g-ix) = 2^+{x), (3.8) 


(3.4) 

(3.5) 

(3.6) 


and using (13.4p . we obtain the identity 
2az) = 2giz)-V{z)-£. 
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(3.9) 









3.2 The RH problem for orthogonal polynomials 


An effective way to characterise orthogonal polynomials with respect to a weight w is 
via the well known RH problem due to Fokas-Its-Kitaev M- In our case the weight 
w{x) takes the form 


w{x) = x" exp 


—n 





(3.10) 


and the relevant RH problem is specified below. 


RH problem for Y 

(a) y : C \ [0, +oo) ^ is analytic. 


(b) The limits of H as z approaches (0, +oo) from above and below exist, are contin¬ 
uous on (0, -|-oo) and are denoted by and IT respectively. Furthermore they 
are related by 

Y+{x) = Y_(x) , X G (0,-koo). (3.11) 


(c) Y(z) = (/ + 0{z as z ^ oo. 

(d) Y{z) = Or^ asz^O. 


This RH problem has a unique solution, 


Y{z) 


f Pn{z) quiz) \ 


(3.12) 


where pj is the degree j monic orthogonal polynomial with respect to the weight w{x), 
hj is the norm of pj , as in (I1.59P , and 


Qj{z) 


J_ pj{x)w{x) 
2'Ki Jq x — z 


(3.13) 


3.3 First transformation Y T 

Define 

T{z) = e-f(3 14) 
The matrix T satisfies the RH problem below. 

RH problem for T 

(a) r : C \ [0, -|-oo) —)• is analytic. 

(b) T has the jump relations 

r+(x) = r_(x)(^ g ^e(o, 6 ), (3.15) 

T+{x) = r_(x) (^Q 1 J > ® e (6,00). (3.16) 
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Outside region 



Yi 

Figure 2: The jump contour in the RH problem for S. 


(c) T{z) = I + 0{z as z —oo. 

(d) T{z) = 0(l)e-2 '"3 as ^ ^ 0 . 


3.4 Second transformation T i— )■ S' 


We now perform the standard technique of opening the lens. This involves introducing 
a contour S5 = (6,00) 7^ as shown in Figure [2j As can be seen in the figure 

the contour divides the complex plane into three disjoint regions. We refer to the 
finite regions in the upper and lower complex half plane as the upper and lower lens 
respectively. We then dehne 


S{z) 


' T{z), 

< 1) ’ 


for outside the lens, 
for 2; in the upper lens, 

for 2; in the lower lens. 


( 3 . 17 ) 


It is then a simple matter to verify that S satisfies the following conditions: 


RH problem for S 

(a) S : C \ S 5 — >■ is analytic. 

(b) S has the jump relations 

5+(z) = S-{z) 

5+(z) = S_(2:)(_°_„ 

/l 2,°'p2ni{z) 

S+{z) = S.{z)[^^ 

(c) S{z) = I + 0{z~^) as 2; ^ 00. 



2; G 71 U 73, 

( 3 . 18 ) 

^ G 72, 

( 3 . 19 ) 

2; G (6, 00). 

( 3 . 20 ) 
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(d) The asymptotic behaviour of S near zero is 

lens, 
lens, 

lens. 
(3.21) 

Remark 3.1 It is important to note that on the contour 71 U 73 U ( 6 , 00 ) away from the 
origin and b, the jump matrices tend to the identity matrix exponentially as n —)• 00 , 
t —>• 0. This suggests we may approximate S' by a global parametrix throughout the 
complex plane minus small neighbourhoods around the origin and b in which it will be 
necessary to construct local parametrices. 

3.5 The global parametrix: 

Following on from Remark 13.11 a first approximation to S can be obtained by ignoring 
the entries of the jump matrices which are exponentially suppressed as n —)• 00 . We 
define the global parametrix (z) as a solntion to the RH problem obtained in this 
manner. 


I, 


Ti / i \ ^ 

S{z) = 0 {l)e- 2 {^) <"3 X 


1 0 
—1 
1 0 


for z outside the 
for in the upper 

for 2 : in the lower 


RH problem for 

(a) : C \ [0, 6 ] ^ is analytic. 

(b) has the jump relations 

= ze(0,6). (3.22) 

(c) = I + 0{z~^) as 2 : ^ 00 . 


Note that this RH problem does not define a solution uniqnely, since nothing is imposed 
on the behaviour of near 0 and b. A solution to the above RH problem can be 
obtained by defining 


P^°°\z) = 2-^^^b^N- 


where N = 2 2 (/ -|- iui), and 


2;2 


D{z) : = 

^{2z-1)2‘ 

in which the branch is chosen 


^ (~ 7 ~) ^ ( 3 - 23 ) 

where D is the Szego function 

(p{z) := z + 'Jz^ — (3.24) 

such that '^{z) ~ 22 as 2 —)• 00 and with cut on [— 1 , 1 ]. 
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3.6 Local parametrix near b 

Given that the only difference between the RH problem for S appearing here and that 
of [26] is the local condition (d) ()3.2ip near the hard edge 0, we see that the local 
parametrix required in a small disk around the soft edge b will satisfy exactly the 
same RH problem (up to a rescaling by b) appearing in [261 Section 3.6]. As is shown in 
[261 Section 3.6], and is by now very familiar, the local parametrix near the soft-edge at 
b can be built using the Airy model problem, and it satisfies the following RH problem. 

RH problem for P^^'> 

(a) \ 'Ps is analytic. 

(b) has the same jumps as S restricted to the contours n S 5 . 

(c) P^^\z) = P^°°\z){I + as n —)• 00 uniformly for z G 

(d) S{z)P^^\z)~^ is analytic at z = 6 . 

3.7 Local parametrix near the origin 

We also require a local approximation to S in a neighbourhood of the origin. This 
will take the form of a local parametrix P^^^ valid within a disk of fixed size, 
characterised by the following RH conditions. 

RH problem for P^^^ 

(a) P^^\z) : \ Ss —^ is analytic. 

(b) has the same jumps as S restricted to the contours 17^°^ n S 5 . 

(c) P^^\z) = + 0 {n~^)) uniformly for z G 9?7®, in the double scaling 

limit where n —)• 00 and simultaneously t —)• 0 in such a way that t —>■ 

s G (0, -t-00). 

(d) As 2 ^ 0 , 


1 


p(o)(2) = 0(l)e"t(l)'''"3 


( ^ 

1 J 

y^-a^-2naz) I J 


for z outside the lens 
for z inside the upper lens 

for z inside the lower lens 

(3.25) 


A convenient way to construct P^^'^ is to first transform the problem to one with con¬ 
stant jump matrices. This can be achieved by the transformation 

P(°)(z) :=p(°)(z)(-z)3"^3e<(^y3. (3.26) 

Applying this transformation, we obtain a RH problem for P^^^. 


23 





RH problem for 

(a) \ Tjs —)• is analytic. This follows from the properties of P^^'^ 

together with the fact that ^+( 2 ^) “ for x G (—oo,0). 

(b) P^^^ has the jump relations, 

pf\z) = p[^\z)(^^^ J), ze 72, (3.27) 

pf\z) = P^\z) (^Ji^ , z e 71 , (3.28) 

Pf (z) = P®(z)(^_L ^G 73. (3.29) 


(c) P^^\z) = P^°°\z){I + 0{n ^))[—z)2^^3Qn^{z)(T3 uniformly for 2 : £ dU^^\ in the 
double scaling limit where n ^ 00 and simultaneously t —>• 0 in such a way that 
2 ^'^cin k t^sG{0,+oo). 

(d) As 2 ; —)• 0, 


p(0)(2) = C)(l)e-i(l)"‘"3(_2)5“'"3 , 



for 2 outside the lens, 
for 2 inside the upper lens, 

for 2 inside the lower lens. 

(3.30) 


We now make the connection to the model problem for T(<^,s) by asserting that 
the function P^^^ can be expressed as 

P^^\z,t) = E{z)^ (n‘^f{z),2~icin~^ 0 - 3 , (3.31) 

where ci = bh{0)^ > 0 and E{z) and f{z) are functions which we will now define. We 
take / to be 

^ i (/ Hs)Ris)ds^ . (3.32) 

Then / is a conformal map from a neighbourhood of 0 to a neighbourhood of 0, with 
/'(O) = —Cl, and we have = (—l)"'e"'^(^). This follows from the definition (13.71) 

of ^ together with the fact that 'tl!{x)dx = 1. Let 

_ 1 2k+l 

Sn,t = 2 kcin k t, (3.33) 

and assume that t depends on n in such a way that Sn,t lies in a compact subset of 
( 0 , + 00 ) for n sufficiently large. 

It can be checked in a straightforward way that ()3.3ip satisfies conditions (a), (b) 
and (d) of the RH problem for Pq. To satisfy property (c), we take 

E{z) = (-1)”P(°°1 ( 2 )^ 73 (nV( 2 ))^ . (3.34) 
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Using the jump relation for and choosing the branch cuts for and 

on the positive half-line, one shows that E(z) is indeed analytic in \{0}, 
with a singularity at 0 which is removable. See [2^ for further details. For z S dU^^\ 
we have that 'n?f{z) is large, and using the asymptotics (jl. 2 ip of 'h at infinity, which 
are uniform for s in compact subsets of ( 0 , -|-oo), one easily obtains condition (c) of the 
RH problem for P. 


3.8 Final transformation: small norm RH problem 

Define 

' S{z)P(°^\z)-^ 

R{z) = < Siz)P^^\z)-\ 

S{z)P^^\z)-\ 

Then R satisfies the following RH problem 


for z e C \ (UW U C/W U S 5 ) , 
for z G 
for z G 


(3.35) 


RH problem for R 

(a) with 

u u (^Ss \ u U [0, 6])) 

is analytic, with the circles around 0 and b oriented clockwise. 


(b) R satisfies the jump relation R+{z) = R-{z)Jr{z) for z G with 


Jr{z) 


'p{0)(^)p(oo)(^)-l, 

^P^°°\z)Js{z)P^°^\z)-^, 


where J 5 denotes the jump matrix for S. 


for 2 : G dU^^\ 
for 2 G dU^^\ 
for 2 G Sr\ 


(3.36) 


(c) As z ^ 00 , we have 

R(z) = / + ^ + 0(z-2), (3.37) 

where i?i is independent of z but depends on n. 

The above RH conditions are straightforward to verify using the properties of S', P^°°\ 
P^^\ and P^^\ For condition (a), one needs to verify in particular that S'(z)P^^)(z)“^ 
is analytic at 0. To see this, one needs to make use of the local behaviour of 'k at the 
origin, see (ll. 22 h or (11.291) . and of the jump relations for S and P^^\ 

From the matching conditions for P^^^ and P^^^ with and the fact that the 

jump matrix for S is exponentially close to the identity on \ as 

n ^ 00 , we observe that 

Jr{z) = I + 0{n~^), z G S/j, n —>■ 00 . (3.38) 

This estimate holds in the double scaling limit where n —)• 00 and simultaneously t —)■ 0 

_1 /U 2 fe+l 

in such a way that 2 ' cin fc t —>■ s G (0, -|-oo). By the general theory for small-norm 

RH problems [TOl Section 7], this implies that 

R{z) = I + 0{n~^), n —>■ 00 , (3.39) 

uniformly for z off the jump contour S/j in the double scaling limit. 
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3.9 Computation of the error 


The estimate (13.391) will be sufficient for the proof of Theorem 11.101 For the proof 
of Theorem \TM we need more precise asymptotics for R. To that end, let us be 
more precise in ()3.38p . In the remaining part of this section, we will specialise to the 
perturbed Laguerre case where V{x) = x. Then, 6=1 and h{z) = 1. The estimates 
below are valid in the double scaling limit where n ^ oo and simultaneously t 0 in 
such a way that ^ s € (0, Too). 

We can write 

Jr{z) = I-\ — Ji{z) + n —>• oo, (3.40) 

n 

where the matrix is non-zero only on and We will express the first 

correction term in the large n expansion for R, given by 

R{z) = I + R^^\z)n~^ + 0{n~'^), z £ C\T,r, n ^ oo, (3-41) 

in terms of R. On dU^^\ we can use (I3.26p . (I3.3ip . and the asymptotics (11.211) for T 
to conclude that 


Ji{z) 


2/(z)V2 


p(oo)(^ 


-1 


p{oo)(. 




(3.42) 


for z G . On the Airy disk, R can also be computed using the explicit construction 
of the Airy parametrix and the asymptotic behaviour of the Airy function for large 
arguments. After a lengthy computation which we do not include here, we obtain 


R{z) 


-P(°°)(z) 

24/(z)3/2 



p(°o)(2)-i, 


(3.43) 


for z G where / is a conformal map sending 1 to 0, with f'{l) > 0. Substituting 

the large n expansion (13.401) in the jump relation = R^Jr, and noting that R{z) —)■ I 
as z ^ oo, we obtain 


RH problem for R^^'> 

(a) R^^'^ : C \ (5t/® U —)■ is analytic, 

(b) R^+\z) = R^yiz) + R{z) for z G U 

(c) R^^\z) —)• 0 as 2 : ^ 00 . 

This scalar additive RH problem for R^^'^ can be solved explicitly, and the solution is 
given by 


= R f 

2m Jqi 


MO 


dO 


ldU(°)udUW i — Z 

where the integral is taken entry-wise. By (13.371) and (I3.4ip . we obtain 


Ri = -^ [ R{C)dC + 0{n-^^ 

2mn JduWuduW 


n —>■ 00 , 


(3.44) 


(3.45) 
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and the term of order n ^ in this expression can be easily computed using residue 
calculus. In the Laguerre case where V{x) = x, we have as n —)• oo, 


= 


n 


(12a^ - 24a + 11) 


^(12a^ +24a + ll) 

_ 1 _|_ ^ 

16 4 . 


2n 


-1 


+0(n-2). (3.46) 


The first term comes from the Airy parametrix, the second term is the contribution 
of the parametrix near 0 , and contains the function r related to the fe-th higher order 
PHI equation. 


4 Asymptotic analysis for 4^ as s —+oo 

In this section we compute the asymptotic behaviour for 4 '( 2 :, s) as s —)■ oo. As a 
corollary, we also obtain the asymptotic behaviour as s —)■ oo of the relevant solution to 
the A;-th PHI equation. Our approach in this section is based on a standard steepest- 
descent analysis. 

4.1 Re-scaling of the model problem 

We begin by performing a transformation T i—)• T which rescales the variable z in a 
convenient way. We dehne 

k 2k 

T{z',s) = z;s). (4.1) 

Then we have the following RH problem for T. 


RH problem for T 

(a) r : C \ S —is analytic, with S the jump contour for T. 

(b) T has the same jump relations as T, 


r+(z) = r_(z)(J 
T+{z) = T.{z)(^_y^ 
r+(z) = r_(z)(_^l„„ J), 

(c) As z —)• oo, r has the asymptotic behaviour 

r(z)= (/ + 


Ti{s) 


+ 0{z 


- 2 ' 




Ne 


S 2fc+l 


Z € 


z € Si, 


•2 G S 3 . 




(4.2) 

(4.3) 

(4.4) 


(4.5) 


where 


fc (j, f qis) —ir{s)\ _ 

Oi’O-l I 21 V / ^ / 1 AlrS.0 


Tl{s) = S 2fc+l5 4fe+2 


ip{s) -q{s) 


g 4fc+2'^3. 


(4.6) 


(d) As 2 ^ 0, there exists a matrix To(s), independent of z, such that T has the 
asymptotic behaviour 


T{z) = Tois)iI + 0{z))e 


(- 1 ) 


(4.7) 


for z G Qj. 
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^3 

Figure 3: The deformed jump contours appearing in the Riemann-Hilbert problem for 
S. The dashed lines denote the location of the original contours Si and S 3 . 


4.2 Contour deformation 

We now note that we can deform the contours Si and S3 to any contour in the up¬ 
per/lower half plane. In particular we deform Si and S3 to the union of {zq, 0) and S'^ 
and S 3 respectively, where zq < 0 is a fixed point which we specify later, and where 
S'^, S 3 are shown in Figure [3j Recall that the angle between S 3 , S 3 (as for Si, S3) and 
the real line can be chosen freely between 0 and vr. 

We denote the collection of contours by S|^ = U S 3 U S 3 . We see from the figure 
that this deformation results in new regions n'-. For z in region for j = 1,2,3, we 
define S to be the analytic continuation of T from region to region fl)-. 

RH problem for S 

(a) S : C \ S'^ —>• is analytic. 

(b) S has the jump relations 

5+(z) = S_(z)(j 

s+{z)=s.{z) j), 

5+(z) = 5_(z)(_^L„ y, 
s+{z) = s.{z) y ^ 

(c) As z ^ 00 , S' has the asymptotic behaviour 

^ + 0(z-^)) (412) 


5(.) = (/+hI 


z G (- 00 , Zo), 

(4.8) 

z G S 3 , 

(4.9) 

z G S 3 , 

(4.10) 

z G (zo,0). 

(4.11) 
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(d) As z —)• 0, 


2k 

/ 1 ^A: + l s 2fc + l _ ^ 

S{z) = To{s){I + 0{z))e^~^^ 


(4.13) 


Note that the only remaining sector around the origin is now and this is the reason 
why Hj in RH condition (d) is replaced by Hi = I. 


4.3 The ^(-function 

We define the g'-function as follows, 

g(z) = {z- zof^'^Pk-i{z)z~’", (4.14) 

where {z — zq)^^"^ is defined to be analytic on C \ (oo,^;o] and positive for z > zq, and 
where pt-i is a monic polynomial of degree k — 1 fixed by the condition 

g{z) = {-l)^+^z-’^ + 0{1), z^O. (4.15) 

This fixes the polynomial pk-i and the constant zq, we have 


k-2 


Pk-i{z) = / ^ + '^PjZ\ 
1=0 


2 ij! ' 


/3, = i-iy+^-\-zo)-^-^ 
zo — 


{2k-l)\\ 

where {2j + 1)!! = (2j + l)(2j — 1)... 3.1. Note that we have 
g{z) = z2 + giz~^^^ + g 2 Z~^^‘^ + 0{z~^''^), z ^ oo, 


where 


gi = h-2 - 2^0, 

3 3 

92 = -^zl + Pk-3 - ^ZoPk-2- 


(4.16) 

(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 


We end this subsection by proving the following proposition which will be needed to 
show that the jump matrices for U (defined below) tend to / as s —)• +oo, for z away 
from the origin. 

Proposition 4.1 (i) For zq < z < 0, we have Re{g{z)) < 0. 

(ii) There exists a 0 G (0, vr) such that for 9 < \ avg{z — zq)! < tT; Re (giz)) > 0. 
Proof. We rewrite (|4.14p as 

(4,22) 
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where 


k-l 

Pk-i{z) := 

j=0 


(2j±l)!!^, 
21 j\ 


(4.23) 


Now note that for 0 < z/zq < 1, g{z) < 0, thereby proving the hrst part of the proposi¬ 
tion. The second part follows from the Cauchy-Riemann conditions if we demonstrate 
that lmg'_^_{z) < 0 for z < zq. To this end we compute 


g'U) 






where 


qk-i{z) 


-k 




i=i 


21 j\ 


(4.24) 


(4.25) 


Due to the negative coefficients in q, we have indeed that lm.g'j^{z) < 0 for z/zq > 1, 
from which the proposition follows. □ 


4.4 Normalisation at oo and 0 of the RH problem 

We now use the ^-function to normalise the behaviour of the RH solution near the 
origin and at inhnity. Define 

U{z-s) = [I+ is^gia-)S{z-s)e-^^9{zy\ (4.26) 

This transformation leaves the jump contours unchanged and gives us the following RH 
problem for U. 


RH problem for U 

(a) U : C \ S'g —)• is analytic. 

(b) U has the jump relations 


t/+(z) = C/_(z)(j 
U+iz) = U.iz) 


1 


0 ^ 


_g7riag-2s 2fc+lg(2) ^ I 

/ 1 0 ^ 


(c) As z ^ oo, U has the asymptotic behaviour 
U{z) = (l+ + C>iz-^)] 


Z e (-00,Zo), 

(4.27) 

z ^ Til, 

(4.28) 

z G S 3 , 

(4.29) 

z G (zo,0). 

(4.30) 


(4.31) 
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with 


Ui{s) =(/ + is^>^+i gia-)x (4.32) 

k 

Ti{s){I — is^>=+^ gia-)+ 

X 2k k 1 3^ 

—s2fc+i g?I + is2fc+ig^(T+ — i(s2fc+TH—s2fc+i gf)a- 
2 6 

(d) As z —)• 0, 

Uiz) = Uo{s)iI + 0{z))z'^^\ (4.33) 

Note that we have 

{Ui)i2{s) = is2^gi{s) + (ri)i2(s), (4.34) 

and hence, by (j4.6ll . 

k 2k 

r{s) = is 2^+1 (C/i)i 2 (s) + s 2 fc+i 5i(s). (4.35) 

By Proposition 14.11 if we choose S 2 close enough to (— 00 , zq), we have that the 
jump matrices for U converge to the identity as s —>■ +00 on the part of the jump 
contour away from the real line. The off-diagonal entry of the jump matrix on (zqjO) 
also tends to 0. The convergence is not uniform near zq, and therefore we need to 
construct a local Airy parametrix near zq to model the jumps near zq, and a global 
parametrix to model the jumps on the real line. 


4.5 Global parametrix 

Ignoring exponentially small jumps and a small neighbourhood of zq, we have the 
following RH problem which we will solve explicitly. 


RH problem for 

(a) : C \ —)• is analytic. 

(b) has the jump relations 

pj°°)( 2 :) = P^°°\z) , z£ (- 00 , 2 : 0 ), (4.36) 

P^\z) = ze{zo,0). (4.37) 

(c) As 2 ; ^ 00 , has the asymptotic behaviour 

p(^){z) = (7 + + 0(z-^)'j z-^^^N. (4.38) 

(d) As z ^ 0, has the asymptotic behaviour 

p(°o)(^) = 0(l)zt‘"T (4_39) 


The solution to this RH problem can be constructed from the one in m and is given 
by 


p(oo)(^) ^ 


1 

ia{-zo)^ 


{z - zq) i N 


_ii3 ( ^z - zo + {-Zq) 2 


^y/z - Zq - (-Zo)2 


(4.40) 


31 






4.6 Local parametrix near zq 

Using the model RH problem for the Airy function, a solution to the following RH 
problem can be constructed in a fixed size disk D around zq. 


RH problem for pUo) 

(a) pRo) is analytic in P \ 

(b) For z & D f] Sg, p(^o) satisfies exactly the same jump relations as U. 

(c) For z S dD, pUo) matches with the global parametrix, 

pUo)(_2) = pM(2)(/ + 0(s"^)), s^+oo. (4.41) 

(d) is analytic at zq. 

The construction of p(^o) proceeds in the standard way. We first introduce the matrix¬ 
valued function 


T{z) 


Ma X < 


( Ai (z) 
VAi'(^) 

( Ai {z) 

^Ai'(^) 

( Ai {z) 

[Ai'iz) 

( Ai (^) 

VAi'(^) 


Aiico^z) \ 
(z^Ai'ico'^z)) 

Ai (uj^z) \ 
io^Ai'iio^z)) 

—w^Ai {ujz)\ 
—Ai'{ujz) j 

—uPAi {oJz)\ 

—Ai'{ujz) J 



\^1 



for z £ 14yip, 

for z £ 

(4.42) 

for 2 : £ flA, 3 , 
for 2 £ flA, 4 , 


with a; = e 3 , Ai the Airy function. 


Ma := (J “.) , (4.43) 

and the regions 14pp dehned as in Figure [H The function T solves the following RH 
problem. 


Airy model RH problem 

(a) T : C \ Sp —>■ is analytic where Sp = 

(b) T has the jump relations 

T+(2) = T_(2)(_^^ j), ^£Sa,i, (4.44) 

T+(z) = T_(z)(j z£Sp,2, (4.45) 

T+(z) = T_(z)(_^^ J), ^GSp,3, (4.46) 

T+(2) = T_(2)(J ze^AA- (4.47) 
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Figure 4: The jump contour T,a for the Airy model problem. 


(c) As 2 : —>■ oo, we have 

T{z) = z-^N (^I + 0 ( 2 - 1 ) 


(4.48) 


Using the Airy model problem we can then write as 
pUo)(^) = ^i(2)T(s^/i(2))e-*^3Ufo30(2), 
where Q{z) = for ilmz > 0 , and /i is given by 

—k 


fli^) ■= n(-^o) M 1 “ 


ZO 


Pk-1 


(4.49) 


(4.50) 


Note that the specific form of fi has been chosen to satisfy —|(s®'“+3 /)2 = 52 ^+ 1 ^. 
Furthermore, using properties of pk-i it is straightforward to show that /i is conformal 
in D when D is sufficiently small, with /{{zq) > 0. Finally, we let 


Ei{z) = P^’^\z)e{z)-^N-\s^fi{z))^. 


(4.51) 


Again it is straightforward to demonstrate that Ei is analytic in D thus showing that 
(I4.49|) possesses the correct jumps. The specific form of E{z) is chosen as to satisfy the 
boundary condition (c) in the RH problem for pUo)_ 


4.7 Small norm RH problem 

Define 


R{z) 


[/(2)p(oo)(^)-i^ forzeC\D, 
[/(^)pUo)(2)-i, foi zGD. 


(4.52) 
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Then R is the solution to a RH problem normalised at infinity and with jump matrices 
uniformly close to the identity as s —)• +oo. One shows that 

R(z,s) = I + 0(s~^), -^R(z,s) = 0(s~^~^), (4.53) 

as 

uniformly for z off the jump contour for R, as s ^ +oo. Defining Ri in terms of the 
large 2 ; expansion of R, 

Riz) = I+^ + Oiz-^), (4.54) 

_ k 

we have in particular that Ri{s) = 0(s ^fc+i), 

4.8 Asymptotics for r(s) and y{s) 

For z large, we have 

U{z) = R{z)P^°°\z). (4.55) 

Comparing the large 2 ; expansions (I4.3ip . (I4.38p . and (14.541) for [/, and R, we 

obtain 

Ui{s) = P^{s)+Ri{s) = P^{s) + Ois~^), s^oo. (4.56) 

It follows from (14.351) . (I4.20p . and the large z expansion of given by (14.401) that 

2/c 3 1 ^ 

r{s) = s2w(/3fc_2 - -zo) - (- 20 ) 252 ^ 0 ; + 0(1), (4.57) 

as s ^ 00 . By (I2.10p . we obtain pi.331) . 

5 Asymptotic analysis for as s —0 

In this section, we compute the asymptotic behaviour of 'I'( 2 ;,s) as s —>■ 0. As a 
corollary, we obtain the asymptotic behaviour as s —)• 0 of the associated solution to 
the A:-th Painleve III equation. 

5.1 The Bessel model RH problem 

If we set s = 0 in the RH problem for T, we obtain a RH problem which is equivalent 
to the Bessel model RH problem used in [Ml EH- This suggests that the solution to 
the Bessel model problem will be a good approximation to T. For 2 ; close to the origin, 
this is not true anymore because of the singular behaviour of T near 0 for s 7 ^ 0 , which 
is not present for s = 0. This is the reason why we will need to construct a local 
parametrix near the origin. 

Bessel RH for T 

(a) T : C \ S —)• analytic, with 

£ = £1 U £2 U £3 
as in Figured] 
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(b) T has the jump relations 
T+(z) = T_(z)(j 
T+(z) = T_(z) 


1 0 

jTrio. 


T+W = T_(2) ( _^ho J 


(c) As z ^ oo, T has the asymptotic behaviour 

T{z) = {l + 0{z-^)) 
where + *< 71 )- 

(d) Conditions as 2 ; ^ 0; 

(i) If a < 0 then T{z) = O 

(ii) If a = 0 then T{z) = O 


z G S 2 , 


z G El, 


2 ; G S 


3- 


(5.1) 


(5.2) 


(5.3) 


(5.4) 


I I — I I — 

\z\ 2 \z\ 2 

\z\ 2 \z\ 2 

log \z\ log I 2 ; 
log| 2 ;| \og\z\)' 


(hi) If a > 0 then 


T(z) = 


O 

O 


\z\ 2 

\z\ 2 

_ ^ 

\z\ 2 

I I _ ^ 

Z 2 


\Z\ 2 \ 

I 1 —S I 5 
\z\ 2 / 


\Z\ 2 

_ ^ 

\z\ 2 


for 2 : G Hi, 
for 2 : G 02 U II 3 . 


(5.5) 


The solution to the above RH problem can be constructed using Bessel functions |26l 
EH, we have 




T{z) = ( / + o(4a"^ + 3)(T_ tt 


niz^r^iz^) -z2K'^{z2) 


H,. 


(5.6) 


when z £ Qj. In the construction of the local parametrix for T near the origin, we will 
also require more detailed knowledge of the behaviour of T at the origin. It can be 
verified that for 2 ; G Ilj and a ^ Z, we have 


T(2:) = To(2;)2;'^ 
whereas for a G Z, we have 

T{z) = To{z)z 2 I ttV j ^ ^ -2 \ Hj. 


(5.7) 


(5.8) 


In both cases To( 2 :) is an entire function. 


35 




5.2 The local parametrix near the origin 

Around the origin, T ceases to be a good approximation to 'h due to the singular 
behaviour of 'h. We therefore construct a local parametrix F around the origin in a 
disk of the form D = {z £ C \z\ < e} with e > 0 hxed but sufficiently small. The 
construction done here deviates from the construction done in m- We prefer this 
modified construction because it is more natural given the behaviour of near 0 given 
in (fL2^ . 


RH problem for F 

(a) F : D \ S —)• is analytic. 

(b) F has the same jumps as T on S. 

(c) On the circle \z\ = e we have the following matching conditions as s —)■ 0, 


Fiz,s) 


T( 2 )(/ + 0(s^)) fora^Z, 
T( 2 )(/ + 0 (s^)), foroGZ. 


(5.9) 


(d) The function TF ^ is bounded near z = 0. 

The solution to this problem takes the following form, 




H. 




(5.10) 


where f 2 {z,s) is given in (|1.28l) . Tq in (I5.7I1 - (I5.8|) . Hj in (ll.23l) - (ll.25l) . and 

Hz) = -hH, 0) + -z", if a ^ Z, (5.11) 

iZSlRTTa 

Hz) =-f2{z,0) +-log^, if a G Z. (5.12) 

TT 2 

It is easily verified that h is an entire function. 

Conditions (a) and (b) of the RH problem are direct, since F is of the form 

Analytic function ^ z^H j, 


just like 'k, see (jl.29l) . and therefore F has precisely the same jump conditions as T. 
To verify condition (c), note that for \z\ = e, a ^ Z, by (15.71) . 


F{z)Tiz)-^ 

= ToW(J + (-1 

= / + O(s^), as s —0. 


22sin7rQ 

1 


To(^) 


-1 


(5.13) 


Condition (d) follows from (jl.29p and ()5.10p . In the case a G Z, the RH conditions are 
checked in a similar way. 
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5.3 Small norm RH problem 

We can now construct a small norm RH problem. Define 


R{z, s) 


"^{z,s)T ^{z), for | 2 ;| > e, 
s)F(z, s)“^, for l^l < e. 


Then i? is a small-norm RH problem and we have 
R{z,s) = I + 0{s’^), s^O, 

and 

—R{z,s) = 0{s^~^), s —)■ 0, 

as 

uniformly for z off the jump contour for R. 


(5.14) 


(5.15) 


(5.16) 


5.4 Asymptotics for r(s) as s ^ 0 

To obtain the initial condition for r, we compute (Te”^ ^^N~^z 4 )^2 as z ^ oo in two 
different ways; firstly using (jl.2ip and secondly using the small s asymptotics obtained 
above. This gives 

r(s) =(l — 4a^)-|-O(s^), r'(s) = C)(s^“^), s —0, s > 0, (5.17) 

8 

valid for any a > —1. By (j2.10jl . (j5.15p . and (|5.16p . we obtain (| 1.341) . 

Remark 5.1 Recall the asymptotic condition (ll.2ip for 'I'( 2 ;,s) as z ^ oo. We used 
before that this condition holds uniformly for s in compact subsets of (0,-|-oo). As a 
consequence of the asymptotic analysis as s ^ 0 done in this section, it follows that this 
expansion is uniform even for s small. This implies that the matching condition between 
the global parametrix and the local parametrix near 0 in the asymptotic analysis for 
the orthogonal polynomials holds also for Sn,t small, and that the error term 0{n~‘^) for 
R in (j3.41l) is uniform also as n —>• oo and t —)■ 0 in such a way that t —)• 0. 

We will use this fact later on. 


6 Asymptotics for the kernel near the origin 

The correlation kernel (|1.60l) can be expressed in terms of the matrix Y as 


Kr 


Im X — y 


( 6 . 1 ) 


To prove Theorem 11.101 we use the large n asymptotics obtained from the steepest 
descent analysis in Section [3l By inverting the sequence of transformations Y T 
S' I—>• R, we may express Y[z) for z in the upper lens region as 


Y[z) = R{z)E{z)'^ (n^f{z),2 '=cin k t 


i7r(a. — l) 

. 2 ^3 


} J) w{z) 2 ^^. (6.2) 
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Inserting this into (16.11) . we obtain 


Kn{x, y; t) = (-1, -l)e (nV(y), ^{y) 

X L{x)^+{n^f{x),Sn,t) (1,-1)^, (6.3) 

where L{z) := R{z)E{z) and Sn,t given by (j3.33p . The boundary value is taken 
according to the oriented contour S in Figure [H which means that we take the limit from 
region II 3 . We now zoom in on the origin by defining x = —Ci^n~^u and y = —c^^n~‘^v, 
with (as before) ci = —/'(O), and consider the asymptotics of Kn{x,y,t) in the double 
scaling limit where n — 00 , t —)• 0 in such a way that Sn,t remains bounded. First we 
note that 


L ^{y)L{x) = I + {u — v)0{n ^), n —)• 00 . 


(6.4) 


Additionally we have 

n^/(x) = u (1 + 0(n“^)) , n — 00 , (6.5) 

and similarly for f{y). If Sn,t s as n ^ 00 , t —^ 0, with s € (0, + 00 ), we have 

{n^fiy), Sn,t) ^+iv, s), T+ {n^fix), Sn,t) s). (6.6) 

Substituting (|6.4p . (|6.5p . and (16.61) into (16.3p . and moreover expressing 4'+ in terms of 
the functions V’l and '(/’2 by (ll. 68 h . we obtain 

,. f _!!: _^ . g \ s)'iIj 2 {v-, s) -'ipijv; s)' 4 > 2 {u; s) 

n^oo civ? "■ \cin2 ’ cire^ ’ C 2 n 2 +i/^y 2m{u — v) ’ 

after a straightforward calculation. This proves Theorem 11.101 

Finally, using the asymptotics of Tp^js) as s ^ 00 and s —>• 0, we can prove 
Theorem 11.121 In the case s ^ 0, we have 

T+(z, s) = {I + 0{s^))T+{z) , (6.8) 


uniformly for z bounded away from zero. Now we can substitute the above expression 
and the expression (15.6p for T in terms of Bessel functions into (jl. 68 j) . Using the 
resulting expression for ipi together with the definition of we obtain, in a similar 

way as in m Section 5.2], the first part of the theorem. 

The limit of the kernel as s —)• 00 is more complicated; one must invert the sequence 
of transformations performed in Section 0] which mapped the RH problem for 4' to the 
Airy model RH problem. The result of this is that, for z close to zq, we have the s —>■ 00 
asymptotics 


4'(z, 


s) =s — is 2 R(s ^z,s)Ei(s ^z,s)x 

for z £ Hi U U Hg, 


T(s3/i(s ^z))0(z) 



for z £ n H 3 , 
for z G fl'i n Q 2 - 


(6.9) 


where y = 2k/{2k + I). Using the above expression in (II. 68 I) we can compute asymp¬ 
totics for i/'j as s —)• 00 in the regions Hg and Q'l H H 3 , and in particular on the negative 
real line. Substituting these asymptotics into the definition of we obtain the 

second part of the theorem. 
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7 Asymptotics for the partition function 


In this section, we will prove Theorem 11.91 by deriving large n asymptotics for the 
partition functions in the perturbed LUE and perturbed GUE. 


7.1 The partition function in the pLUE 

The main ingredient required to obtain asymptotics for is a differential identity 

with respect to the perturbation parameter t. More general identities of the form below 
were obtained in but for completeness we give a proof of the differential identity 

relevant to us. 

We note that the function zTr (Y~^ {z)Y'(z)) is not analytic at infinity: indeed it 
has a discontinuity on the real line. However, since the discontinuity becomes exponen¬ 
tially small as 2 : —)• ± 00 , zTr {Y~^{z)Y'{z)) does admit a full asymptotic expansion in 
negative powers of z as z ^ 00 . We have 

2 Tr {Y-\z)Y\z)) = ^ + 0(2-2), 

and denote Res^=oo (^Tr (Y~^{z)Y'{z)^^ := —ci, although strictly speaking this is not 
the residue at infinity, because infinity is not an isolated singularity. The integral 

j, zlY{Y-^{z)Y'{z)a^)dz 

over a large counter-clockwise oriented circle {z : \z\ = M} is not exactly equal to 
27rici, but it tends to 27rici as the radius M tends to infinity. This observation will be 
important in the proof of the following lemma. 


Lemma 7.1 (cf. [3114]') We have 

^ ^Res^=oo {zTt {Y-\z)Y'{z)a3)) . (7.2) 

Proof. Recall (jl.Op . which, after the substitution Xi = becomes 


ypLUE I'n?-\-OLn ypLUE 

^n,k,a ^ ^n.k.a ’ 


n^k,a 


where 


^pLUE 


1 


n'- Jio,+ooY j=i 


(7.3) 


(7.4) 


To evaluate Z^^^ ^ we introduce monic orthogonal polynomials pj satisfying 


poo ^ 

/ Pi{x)pm{x)w{x)dx = hi6im, 'u;(x) = x“ exp 
Jo 


—n [tx + X 


-k 


(7.5) 


Recalling the orthogonality conditions for pn given in (|1.59p . we easily obtain 

Pn{x) = t~'^Pn{tx), hn = (7.6) 

A similar calculation for the Cauchy transform appearing in (|3.12p shows that 

Pn{x)w{x) 


Quiz) = 


1 r 

2vri Jo 


-dx = t 


X — z 


"'qnitz). 


(7.7) 
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Summarising, if we define 


Yiz) 


Pn{z) 


hn—l 


Qn{z) 


Pn-l{z) -J^qn-l{z) 


hn—1 


(7.8) 


we obtain the relationship 

Y{z) = (7.9) 

Expressing as Z^^^ = nr=^o^ taking the logarithmic derivative of both 

sides, we get 


d , TZpLUE 



dw{x) 

~dr~ 


dx, 


(7.10) 


where we have made use of the orthogonality of the polynomials. This expression can 
be written as 


^ log = “^ J xKn{x,x)dx, 

where 


Kn{x,y) 


n—1 

y/w{x)w{y) 

j=0 


Pj{x)pj{y) 


(7.11) 


(7.12) 


The expression for the kernel Kn in terms of Y, similar to (16.11) . together with the jump 
condition for Y implies the identity 

Kn{x,x) = -^ (Tr {Y-\x)Y[{x)a^) - Tr (y_-n 3 ;)ll(x)a 3 )) , (7.13) 

for X G M, which combined with the fact that Tr{Y~^{x)Y'{x)a 3 ) = tTr{Y~^ {tx)Y'{tx)a 3 ), 
gives 

rJ ^ r? 

=^tl "" '(^)n(^)^ 3 ) - Tr(y_-n^)ll(^)<^ 3 )] dx. (7.14) 

To prove the lemma, by (|7.3p . we need to show that the above is equal to 
^Res^=oo (^Tr (y"^( 2 :)y'( 2 ) 0 - 3 )) • 

To see this, consider the integral of zTi (Y~^ {z)Y'{z)as) over the contour C depicted in 
Figure [5j One one hand, this integral is equal to zero because the integrand is analytic 
inside the integration contour. On the other hand, condition (d) of the T-RH problem 
implies that for sufficiently small e, Tr (y“^( 2 )y'( 2 )o- 3 ) is bounded on the circle C 3 , 
uniformly in e. Furthermore, condition (c) of the T-RH problem implies 

Tr {Y-^{z)Y'{z)a 3 ) = 0(z-^), z ^ 00 . 

Therefore, letting e —)• 0, letting the radius of the large circle Ci tend to infinity, and 
letting ( 72 , 04 approach the positive half-line, we arrive at the identity 


- 27rzRes2=oo (^ Tr (y ^(z)Y'(z)cr 3 )) 

roo PCXD 

+ / xTr (y_^^(x)yKx)(T 3 ) dx — / xTr (y_r^(x)y^(x)(T 3 ) dx 

Jo Jo 

Substituting this into (I7.14p . we obtain the result. 


0. (7.15) 
□ 
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Figure 5: The contour of integration C = Ci + C 2 + + C 4 

We may now give the proof of the first part of Theorem 11.91 We proceed by computing 
the residue at infinity in (17.11) . The expression for the matrix Y that is valid away from 
the support of the equilibrium measure (i.e. outside the disks near 0 and b, and outside 
the lens-shaped region) is 

Y(z) = (7_lg) 

To obtain an expression for the residue, we compute the expansion as z —)• 00 of the 
above expression and substitute them into Tx(Y~^Y 'giving as z —>■ 00 , 

Yi{Y-^Y'a^){z) 

= i + 2n(l + Tr(Picr 3 ) - - Tr(Pio- 3 ) -h \ . 

z \ z'^ z z z J 

Here, the quantity /ii arises from the expansion of the ^-function and is the first moment 
of the equilibrium measure, 

:= f xil){x)dx, (7-17) 

Jo 

and the quantities Pi and Pi are defined as the functions of t appearing in the z —>■ 00 
asymptotics of p(°°) and P, 

p(°^)(z) = I + ^+0(z-^), (7.18) 

R(z) = I+:^ + 0(z-^). (7.19) 

z 

In the Laguerre case, we have ii\ = 1, Tr(Picj 3 ) = —2a;, and Pi is given by (I3.46|) . 
Recall from Remark 15.II that the asymptotics for Pi are uniform also when t —^ 0 very 
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rapidly. Substituting all this in the differential identity (|7.2|) . and integrating between 
0 and t, we therefore obtain 


log 


ypLUE 

^n,k 

rypLUE 

^n,k 


( 0 ) 



4a^) — r(2 fen fe*^ + ^), (7.20) 


in the double scaling limit where n —oo, t —)• 0 in such a way that st^n = 2 fen fe t —>■ 

2fc+l _ 

s e (0,+oo). After the change variables ^ = 2 fen fe t, we obtain (ll.54jl . 


7.2 The partition function in the pGUE 

Proposition 7.2 Let be the partition function in the pLUE, defined by (jl.52l) . 

and let he the partition funetion in the pGUE, defined by (ll.SOp . Then we have 

the following identities: 


log^2n,fe.a 

ypGUE 




= (^ + a + ^) 


n log 


n,/c,a+| 

2 n 2n + 2 

+ (n + 1) log 


(7.21) 




fc+i 

'2n +fe 


2 n + 2 


2 n + 1 


^ n,K,Q:+^ 


+ 


2 n + 1 


fc+i 

' 2n + 1 \ ~ 


2 n 


(7.22) 


Proof. Let Pj{x; n, k, a), j = 0,1, ..., be the monic perturbed Laguerre polynomials 
defined by the orthogonality relations 



Pj{x)pi{x)x°'e 



dx 


hj{n, a, k)5je. 


Then it is well-known that 


(7.23) 


77—1 

log ^nlk!a = hj (n, k, a). (7.24) 

j=0 

A similar identity holds for the monic perturbed Hermite polynomials pj(n,k,a), 
j = 0,1,... defined by 



hj{n, k, a)Sji. 


We have 


(7.25) 


77—1 

= ^loghj{n,k,a). (7.26) 

j=0 

Using (I7.23P and p7.25l) and the uniqueness of the orthogonal polynomials defined 
by those orthogonality conditions, it is straightforward to derive the following identities 
(similar identities can be found in [H Appendix B] and |17l Exercise 4.4]): 

P 2 j{u-, 2n, k, a) = pj{u^;n, k,a — 1/2), 

P 2 j+i{u-, 2n, k, a) = upj{u^‘, n,k,a + 1/2). 
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(7.27) 

(7.28) 



















This implies that 


h2j{2n,k,a) = hj{n,k,a - 1-), h2j+i{2n,k,a) = hj{n,k,a + ^). (7.29) 


Using those relations together with (I7.26P and (I7.24p . we obtain 

n—1 n—1 


y]iog/i 2 j (2 

n,k,a) + log h2j+i ( 2 n, k, a) 

(7.30) 

1=0 

1=0 


n—1 

y]log/ij(n, 

1 1 

k,oi- -) + '^\oghj{n,k,a + -) 

(7.31) 

1=0 

1=0 


logz^y^" 1 

n,fc,Q:— 

TlogZ^y^i. 

n,k,a+-^ 

(7.32) 


Similarly, 


n—1 




y^log/i 2 j( 2 n + 1, k, a) + log/i 2 j+i( 2 n + 1, k, 


a) 


j=0 


j=0 


'^loghj{n + -,k,a - -) + '^loghj{n + -,k,a + -) 


j=0 


j=o 


= logA{n + l,n + ^,k,a - ^) + log A{n,n + k, a + ^), 


where 


1 


A{n,N,k,a) = — I A( 3:)2 17 

n'- j[o,+oo)" 


^.j+i A: 


dxj. 


We have 

A{n,N,k,a) = 


jl \ n^-\-na 


^pLUE 

‘^n,k,a 


^ fe+i 

n 


Substituting this into (I7.33P completes the proof. 


(7.33) 


(7.34) 


(7.35) 


□ 


The proof of (I1.55P can now easily be completed by substituting (11.541) into (17.211) - 

(TM) . 
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